FAU-TP3-04/02 
hep-th/0402014 
February 1, 2008 



Phase Diagram of the Gross-Neveu Model: Exact 
Results and Condensed Matter Precursors 



Oliver Schnetz, Michael Thies and Konrad Urlichs^ 

Institut fiir Theoretische Physik III 
Universitat Erlangen-Niirnberg, Erlangen, Germany 



Abstract 

Recently the revised phase diagram of the (large N) Gross-Neveu model in 1+1 di- 
mensions with discrete chiral symmetry has been determined numerically. It features 
three phases, a massless and a massive Fermi gas and a kink-antikink crystal. Here 
we investigate the phase diagram by analytical means, mapping the Dirac-Hartree- 
Fock equation onto the non-relativistic Schrodinger equation with the (single gap) 
Lame potential. It is pointed out that mathematically identical phase diagrams 
appeared in the condensed matter literature some time ago in the context of the 
Peierls-Frohlich model and ferromagnetic superconductors. 
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1 Introduction 



The Gross-Neveu (GN) model in 1+1 dimensions is probably the simplest interacting 
fermionic field theory one can write down, 



(here, i = l...iV is a fiavor index, and the model is defined through the 't Hooft limit — > 
oo, Ng"^ = const. j2]). Yet, as evidenced by more than 900 citations in the hep-archive to 
date, this toy model has turned out to be quite valuable for a variety of physics questions 
(for a recent review and a discussion of subtleties associated with low dimensions, see 
Ref. If one solves it via semi-classical methods following the original 1974 paper, a 
surprising number of phenomena of interest to strong interaction physics unfold. Already 
in Ref. asymptotic freedom, spontaneous breaking of the discrete chiral symmetry 
~^ 75"^) dynamical fermion mass generation and a scalar qq bound state (the cr-meson) 
were demonstrated. The version with continuous chiral symmetry, the two-dimensional 
Nambu-Jona-Lasinio model j3] (NJL2) possesses an additional massless vr- meson. Shortly 
afterwards, massive kink [S] and kink-antikink type baryons were derived analytically 
in the discrete chiral model, whereas the NJL2 model features massless baryons [3 |Hj. 
Promptly it was shown that chiral symmetry gets restored at finite temperature in a 
second order phase transition ^ llOj. 

In contrast to these early works, it took a rather long and winding road to determine the 
full phase structure of the GN model at finite chemical potential and temperature. In 1985, 
the first phase diagram was proposed IJJJ , followed by a number of works elaborating on it 
(see e.g. [121 US])- Its prominent features were two phases (massive and massless quarks), 
separated by a line of first and second order transitions meeting at a tricritical point. In 
a density-temperature plot, a mixed phase would appear in the region of low density and 
temperature. It can be pictured as droplets of chirally restored vacuum containing extra 
quarks, embedded in the symmetry broken vacuum — reminiscent of the MIT bag model 
At the first order phase transition, the droplets fill all space 

This scenario as well as the underlying phase diagram were believed to hold for both 
variants of the GN model (with discrete and continuous chiral symmetry) . Only in 
2000 it was noticed that the phase diagram was hard to reconcile with the known baryon 
spectrum. In the case of the NJL2 model, a thermodynamically more stable solution of the 
mean field equations turned out to be one where the chiral condensate assumes a helical 
form, namely a circle in the {ipip, ipi'y^ip) plane superimposed with a uniform translation 
along the x axis (a "chiral spiral" jIH]). Since the winding number is equal to baryon 
number, one can regard this structure as a caricature of a Skyrme crystal ^^1; a point of 
view emphasized in Ref. |15^. This helical order parameter was subsequently confirmed 
in Ref. [T7j , the (substantially) modified phase diagram of the N JL2 model is discussed in 



As the bosonization technique based on the massless boson is unavailable in the GN 
model with discrete chiral symmetry, it took somewhat longer to construct a satisfactory 
phase diagram for this model. Following a variational calculation which clearly showed 
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Ref. 0. 
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the existence of a crystal ground state at T = and any finite density [TH], the full 
phase diagram was obtained in Ref. jTHI, based on a numerical solution of the Dirac- 
Hartree-Fock equations. The main new feature is the appearance of a third phase, a 
kink-antikink crystal, in addition to the two previously known homogeneous phases. It 
supersedes the mixed phase in the old phase diagram. All transitions are second order, and 
the tricritical point gets transmuted into a Lifshitz point characteristic for the transition 
from an inhomogeneous to a homogeneous ordered phase in condensed matter physics 
fK)\ . The mechanism which drives the spontaneous breakdown of translational invariance 
at finite density was identified as Overhauser effect i21j with gap formation at the Fermi 
surface. More recently, guided by these results, the analytic form of the self-consistent 
scalar potential at T = has been guessed correctly and the self-consistency of the crystal 
ground state could be established analytically [22] ■ Many of the previous numerical results 
can now be written down in a concise, analytical way. Instrumental for this solution was 
the fact that by using an appropriate scalar potential expressed in terms of Jacobi elliptic 
functions, the Dirac-Hartree-Fock equation could be mapped onto a Schrodinger equation 
with the Lame potential which was solved analytically long time ago |23j . 

In the present paper, we extend the work of Ref. [22] to finite temperature. It turns out 
that the ansatz for the self-consistent potential used in [22] is general enough to encompass 
the self-consistent potential at any temperature and chemical potential. Apart from the 
obvious benefit of analytic insight into the phase structure of the GN model, our results 
have proven helpful in uncovering some striking relationship between the GN model and 
certain problems in condensed matter physics, notably in the context of superconductivity. 

This paper is organized as follows: In Sect. |2l we present our two-parameter ansatz for 
the scalar potential and the solution of the Hartree-Fock-Dirac equation. In Sect. [HI we 
minimize the grand canonical potential and show self-consistency at finite temperature 
and chemical potential. We also discuss how to compute thermodynamic observables. In 
Sect, m we map out the phase boundaries by appropriate series expansions, study the 
T = limit and derive an effective action of Ginzburg-Landau type near the multicritical 
point. Finally in Sect.|Hlwe point out the close relationship between the phase diagram 
of the GN model and quasi-one- dimensional condensed matter problems. 

2 Ansatz for the scalar potential and solution of the 
Dirac-Hartree-Fock equation 

We apply the relativistic Hartree-Fock approximation which becomes exact in the large- 
limit (see [3] for details). The starting point for our discussion is the Dirac-Hartree-Fock 
equation. 





with the following choice of 7-matrices, 



7° = -^1 , 



7 



1 



m, 7 
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7°7^ = -0-2 . 



(3) 
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In the Gross-Neveu model with discrete chiral symmetry, S{x) is reaL In terms of the 
upper and lower spinor components (j)± the Dirac equation consists of two coupled equa- 
tions 

which can be decoupled by squaring, 

f dS 



Eqs. (j21ISl) fall into the pattern of supersymmetric (SUSY) quantum mechanics. As ansatz 
for 5* we choose the superpotential of the Lame potential 

S{x) = Ak — ^ , (6) 



where sn, cn = a/1 — sn^, and dn = a/I — K^sn^ are Jacobi elliptic functions to the 
modulus K ^25|. We introduce dimensionless variables 

e = Ax, SiO = Six)/A, uj = E/A (7) 

and find that S has period 2K, where K = K[k] is the complete elliptic integral of the 
first kind. Moreover the scalar potential satisfies the symmetry condition 

~S{i + K) = -~S{i) . (8) 

With the ansatz the second order equation (jSJ becomes 

(-^ + 2'^'sn2(e' '^)) = + '^')0+(O (9) 

the simplest case of the Lame equation (with a single gap). The solutions are well known 

123 EE!, 

0.(O=Ar^i|±^e-H^ (10) 

where H, and Z are the Jacobi eta, theta and zeta function, respectively P^. (Re- 
gretfully there exist different conventions for the argument of the Jacobi zeta function. 
We choose the convention with Z(K) = 0). The parameter a is related to the (reduced) 
energy uo and Bloch momentum p by 

TT 

|^|=dn(a,K), p=-iL{a) + — . (11) 

Note that a and —a supply the two independent solutions for a given energy eigenvalue 
uj. The condition for real valued energy and momentum leads to two energy bands uo'^ = 
... (1 — K^), p = . . . 7r/2K and u"^ = 1 . . . oo, p = 7r/2K . . . oo. The lower band is 
parametrized by a = K + i(K' — rj), the upper band by a = irj, where in both cases 
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rj = . . .K.' = K[Vl — K^]. From Eq. pij) we calculate (E = E[k] is the complete elliptic 
integral of the second kind) 

^ = ± - (^2) 

^(cj2 - 1 + ^:2) (t^2_i) ' 

where the plus sign refers to the upper band, the minus sign to the lower band. We chose 
the sign of a in such a way that the slope of the dispersion relation is positive in both 
bands. The eigenfunction corresponding to —a has quasi-momentum vr/K — p. The full 
dispersion relation is plotted in Fig. [T] 

With 0+ given, follows from the Dirac equation (jU, 

More explicitly, from Eqs. Q and (jH)) we can read off that a simultaneous translation ^ 
^ + K and a discrete chiral transformation ip 'y^ip leaves the Dirac equation invariant. 
Thus -f^ipi^ + K) oc or equivalently c0+(O = + K), c0_(O = + K) for 

some constant c. By iteration we find 

cV+(0 = + 2K) = -e^^^^0+(O . (14) 

Comparison with ()10j) determines c up to a sign ambiguity which can be resolved by an 
explicit calculation, 

040 = -sgn(.;)e^^("V+(e + K) . (15) 

We still have to determine the normalization factor A/" in Eq. (fTn|) . In a continuum 
normalization, the spatially averaged fermion density should be normalized to 1, 

1 /■2K 

Inserting (fTUj) and using an addition theorem for H (note that }l{x+y)il{x—y)/Q'^{x)Q'^{y) 
is elliptic in x and y and thus determined up to an additive constant by its pole structure) 



H(a; + ?/)H(a; — y) tt 



leads to 



= {dn\y,K)-dn\x,K)) (17) 



W .J^y^ o dn^(a,«:)-dn^(e,«:) 



The spatial average is easily computed. We find 



2 'Knyl — T? 



7re2(«)|dn2(«,fi:)-E/K| ■ ^^^^ 



With Eq. (fTT|) we can simplify the expression for to 

la^^-dn^(e^ 



(20) 



Let us now calculate ipip and ip'^ip., the scalar and baryon densities for a single orbit. We 
find 

= + 010+) = -{d^ + 25)|0+p (21) 

and with Eq. ^ 

Notice that every orbit has the ^-dependence of the Lame superpotential. On the other 
hand Eq. (US} implies = + + K)|2. Using Eq. ^ and the addition 

theorem for dn(,^ + K, k) leads to 

iFiom. the normalization condition (ip'^ip) = 1 we deduce the spatial average 

= {S^) = 2- -2E/K. (24) 

This completes the preparations needed for calculating the grand canonical potential 
density. 



3 The grand canonical potential 

In order to determine the phase diagram of the GN model, one has to minimize the 
grand canonical potential in the space of scalar potentials S{x). We will show that the 
two-parameter ansatz © provides a self-consistent solution so that it is sufficient to 
minimize with respect to A and k. In the relativistic Hartree-Fock approximation, the 
grand canonical potential density per flavor is given by (see jH] and references therein) 



1 /-A/a 
^ = -— / dg In 

pvr JO 



(1 + 



1 + e 



+ 



1 



2Ng^X Jo 



dxS (x) 



(25) 



Here, A/2 is the ultra-violet cutoff, A = 2'K/A the spatial period of S{x), and q = Ap. We 
change to lo as integration variable, spelling out the contributions from the two energy 
bands. With Eqs. ^ and (j21D we find 



A 



did + Jduj In 
V 1 



1 + e 



'l3{Au-ii) 



+ 



2iV^ 



(26) 



^From Eq. (fTT|) we extract the behavior of u for large p and find that = a;(g = A/2) = 
K/2A + As^ / K + 0{hr^). Due to the quadratic divergence, we have to keep the 1/A term. 

In the following it is convenient to combine the integral over both energy bands as well 
as over positive and negative energy modes into the real part of a single integral in the 
complex plane. The path of integration is shifted slightly away from the real axis. We find 
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that dp/do; defines a double cover of tlie complex plane and choose the sheet by requiring 
dp/duj > for u; > 1. The gaps and the sign-change in Eq. (jl2|) match the behavior of 
dp/du in the complex plane. Setting 

a = (3A, z/ = , (27) 

we obtain 

T,(3^^ = ^ _ a hm Re / doo , ~ ^ In f 1 + e-"-+^) . (28) 



The physical value of \1/ for given /3 and /i is determined by the minimum of \l/ in the 
two-dimensional space {k,,A) G (0 . . . 1, . . . oo). 

It turns out to be convenient to minimize \1/ in the directions given by two specific linear 
combinations of d/dA and d/dn, corresponding to ad/da + K.d/dK. and ad/da — sd/ds. 
Consequently we introduce 

= a^(l - E/K) [^d-a-l- l^-E/K^'L ' "^'^ 






1-k2-E/K/ a (2-/t2-2E/K)(l ; ^ , 
= .^(l-E/K) (l-K^-W "^1-'^ * 

where we have converted d/dK. and d/ds into the S/Sk terms. The calculation of Fi and 
F2 is straightforward (replace by tu/a in \1/ before differentiation), 

oo+ie 

F2 = -limReydcj|^ ^ = I -ln(l +e-"^+'^) . (30) 

—oo+ie 

Note that the integral in Fi is logarithmically divergent as A ^ 00, whereas the integral 
in F2 is convergent. Hence F2 is a function of a, u, k only. In the range < k < 1, < A, 
minimizing \l/ is equivalent to Fi = F2 = 0. This defines the crystal phase. It exists only 
in a certain region in the (T, /i)-plane (see |19j and Sect. |3)). In this region the value of 
\1/ has to be compared with the (mathematically) co-existing value of \1/ for the massless 
Fermi gas (corresponding to k = or A = 0) and the value of \1/ for the massive Fermi 
gas (corresponding to k = 1). It has been checked numerically that the crystal solution 
(whenever it exists) is thermodynamically most stable. This result is confirmed by an 
analytical investigation near the phase boundaries as well as at T = 0, cf. Sect. |3 

One advantage of our particular choice of direction in which we minimized the potential 
[cf. Eq. (|29|) ] is the fact that Fi = is equivalent to the self-consistency of our ansatz. To 
show this we manipulate the thermal expectation value (note the sign change for negative 
u in Eq. ^ ) 

A/2 

^^^^^^ = IJ^^^^ i,iE-l, + 1 - e/^(-^-.) + l ) ^31) 
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in the way we reformulated ^. With Eqs. (fT^ . ()22|) we obtain 



00+l£ 



Up^)^^ = -Six) hm Re / Auj , ^ :- ^ . (32) 

^ ' TT ^'e^O J ^(^2_l+^2)(^2_l) (e" + 1) ^ ^ 



Applying a partial integration in Eq. (jHOJ we can use Fi to express {ip4')th as follows, 



(^^)th = ls{x) - ^ ) • (33) 



This in turn reduces to the self-consistency condition S{x) = —Ng'^{iptp)iYi for Fi = 0. 

In order to isolate the divergencies and to make expansions for k ^ 0, k ^ 1, and a — 
easier accessible, we cast \E',Fi,F2 into a different form. The transformation amounts to 
changing from the integrals over the two energy bands to an integration over the single 
gap. Here we simply give the results which can be verified by comparing the Taylor 
expansions of corresponding functions at a = 0. Whereas the results for the transformed 
expressions are trivial (cf. Eqs. (jHH), (j^ )- some work is needed to expand Eqs. (j^ . (jHUj) 
at a = 0. The calculations are lengthy but straightforward if one proceeds as follows: 

1. replace uj hj uj/a as integration variable. 

2. show that the expansion in a = can be derived from expanding the integrand. 



UJ~U 

2 



-2 



3. use integration by parts twice to obtain the factor [ln(l +e '^+'^)]" = [2 cosh 
This renders the integrals finite for A — > cxo. 

4. The a°- term can be converted into a sum and thus evaluated. 

5. Using the residue theorem, the remaining terms can be converted into a sum over 
Matsubara frequencies (the a? term needs some extra attention). 

6. The sum over the Matsubara frequencies can be evaluated yielding rational multiples 
of C(2/c + 1) (where C, is the Riemann zeta function and /c = 1, 2, . . .). 

7. Comparison with the power series of ImlnF Q + at a = allows one to give 
the result in terms of Im In F and its derivatives. 

Proving the results is now equivalent to deriving Eq. 

We reintroduce A and ii to show the dependence on the basic parameters of our ansatz: 

_A2 _ M _ ^ _ Td ^A^ - 

Stt 27r 6/32 27r ^ 7r2 y I ^ K J 







vr , r(l + f (/i + AA^)) A/3 ^ 



'i/5AA^ r(i + f (/i-AA^)) 47r Ng^ 



(34) 



where 

\ = ^ll- K^sin\ip). (35) 

The quadratically and hnearly divergent terms (energy density — A^/Svr and baryon den- 
sity — /iA/27r of the Dirac sea) are irrelevant and can simply be dropped. From now on, 
we will work in units where the dynamical fermion mass in the vacuum is 1. In these 
units the vacuum gap equation reads 

In A. (36) 



Ng 



Invoking this equation, all remaining infinities are removed and we obtain the renormalized 
expression 

6 2 2 te / V ' KaJ r(l + i(i/-aA 



Similarly we obtain {if) = (In F)') 



7r/2 
(3 I f 

-Fl,ren = - lu — + - / dy? Rc 



(37) 



(38) 



da a 



We immediately see that \l/ren(/9, o^, z^^), Fi^^^niP, J^^, '^^) and -^2(0^, z/^, k^) are even, 
analytical functions in a, v and k. A practical procedure to solve -Fi,rcn = -F2 = is as 
follows: 



1. Choose values for two out of the three parameters a = 0...oo, z/ = 0...oo, k = 
0...1. 

2. Solve F2 = (Eq. ()39p ) for the third parameter. This requires finding the zeros of 
a function of one variable given as a one-dimensional integral. 

3. Solve (trivially) Fi^^cn = (Eq. ^) for /3. 

4. With {A = a/jS^n = p/(3,h,T = given, the potential density (jT7|) and all 
observables are given at one point in the crystal phase in the (T, /x)-plane. 



Note that the translationally invariant solution is obtained by setting k = 1 and skipping 
step 2. 

This enables us to compute the phase diagram and to explore all limits of interest (see 
Sect. 13]). We illustrate the result of the minimization procedure by plotting lines of 
constant A and k in the (/i, T)-plane, see Fig. |21 To clarify this diagram, we also show 
separately the lines A = const, and k = const, in Figs. 01 and ^ The thick lines in Fig. |21 



9 



correspond to A = (phase boundary between massive and massless phases) and A = 1. 
Like a "separatrix" the latter curve divides the plot into regions where the lines A = 
const, emanate from the /x-axis (for A > 1) and from the T-axis (for A < 1), respectively. 
The other end of the curves lies on the phase boundary between crystal and massless 
phases. For A < 1, the kink visible in the contour lines reflects the phase boundary 
between crystal and massive phases. In Fig. HI the thick curves marked k = and n = 1 
correpond to phase boundaries between crystal phase and the two others. We shall return 
to this phase diagram in Sect. El 

For later use notice that the potential density at its minimum is equivalent to the simpler 
expression 

7r/3'^i = 7T(3^^,,n-^F, + a''F2 (40) 



2 

7r/2 

dip cos(2<y9)Re 



27r 







The standard thermodynamic observables can be computed in a straightforward manner, 
once the grand canonical potential density (see Eq. ()37|l ) is known. The pressure P, 
the (spatially averaged) baryon density p, the entropy density s and the energy density u 
are given by 

s = ;3^^*, u = Ts-P + i^p. (41) 

The x-dependence of the baryon density in the crystal demands some extra attention. In 
analogy to Eqs. (jSH), we find for p(x) = ('?/^V)th 

A/2 

If.,.,/ 1 1 



p{x) = — J dpip'^ip ^ 





g/3(a;-/i) _|_ ]_ g/3(-w-^t) _|_ ^ 



co+ie 



-limRe / du , ' -?74 ^ 7 42 

TT.-O J ^(^2_l + ^2)(^2_l)e/^(^-A') + l 



A/3 



With Eqs. pHl and (j37jl we readily check that its spatial average satisfies {p{x)) 
—d'^ /dp. We will evaluate p{x) analytically at T = in the following section. 
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4 Phase boundaries, zero temperature limit and mul- 
ticritical point 



Limit K ^ 0, perturbative phase boundary 

Here we are interested in tlie boundary between tlie cliirally restored pliase and tlie crystal. 
In Ref. it has already been determined via almost degenerate perturbation theory. 
As a cross-check we rederive it from the full thermodynamic potential. A straightforward 
expansion of F2 [use Eq. (j^Hl)] yields 



-Re 



+ 0{^) . (43) 



32 da 

With Eq. (jHEI) we can summarize the equations for the k — > phase boundary by 



In = - min Re 

AttT 2 a>0 



+ + «)) + (« 



(44) 



The resulting curve is included in Figs. El and IHl (label "k = 0"). For small u where 
Re ijj (2,^ + < {ip{2, x) = d'^'i/j(x)/dx'^ etc.), the unique minimum is at a = 0. In 
this range of u the phase boundary does not touch the crystal region in the (/x, T)-diagram. 
It corresponds to the transition between the massless and the massive homogeneous so- 
lutions described by 

\n^ = Re^(- + —) . (45) 



4nT \2 27iJ 

The critical temperature at /x = (using ip{l/2) = — C — In 4) coincides with the known 
value 

e° 

T, = - . (46) 

71 

The tricritical point is situated at z/j with 

Re^ (2,1+'^) =0, z/t = 1.910668 (47) 
\ 2 271 J 

One finds y^t = 0.608221, Tt = 0.318329, A = 3.141401 in agreement with the tricritical 
point of the old phase diagram [TT] . 

For V > Vt the right hand side of Eq. ()44|1 develops a minimum for a > 0, whereas a = 
becomes a maximum (implying T(a > 0) > T(a = 0)). Here Eq. (jl^ defines the boundary 
between the massless homogeneous and the crystal phase. In the neighbourhood of this 
boundary a leading order expansion allows us to analytically compare the values of the 
minimized potential density in the crystal phase with the massless homogeneous potential 
\E'(k = 0) = — 7r/6/9^ — /i^/27r. To this end we fix a z/ = y/T = const, line in the {fi,T)- 
diagram. This line has a unique intersection point Pq = (/xq, ^0) with the phase boundary. 
The value oq of a at Pq is determined by the location of the minimum in Eq. ()44|) . Now 
we determine the minimum value of ^'rcn on the u = const, line in the neighbourhood of 
Po as a function of T (and u, qq, Tq) and compare it with the k = result. We obtain 

^ _^ , 16vr(r - To)' 

'±'ren,min — ^vcn,K=0 H r 7 ; — -\ ; "T . [^o) 



Re 



V' (2, 1 + ^{J^ + ao)) + (cto -cto) 
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The minimum condition in Eq. ()44|) implies < d^(r.h.s) / da^ at a = Oq and thus 
^ren,K=o T < Tq, u > Ut. The Crystal solution is thermodynamically favorable. This 
proves that a phase transition occurs. The (T — To)-dependence shows that it is of second 
order. 

Limit K ^ 1, non-perturbative phase boundary 

This limit is relevant for the phase transition between crystal and massive phases as well 
as for the translationally unbroken, chirally broken "massive" phase alone. 

The relation between a and v along the (non-perturbative) phase boundary is determined 
by F2 = or equivalently 

1 / \ (9 1 ''/■^ 1 r + 77-(zz + acoso^)) 

— limlnfl-zt^WafK) = ^- / d(^^— Imln — >^ — — f = . (49) 

4a«-i V ) ^' daa J ^cos(^ r (i + f(z/ - a cos (/^)) 



V2 27r 

With these values for a and u 

F^{k = 1) = -In^ + - r difiRetp (- + —(u + acosip)] = (50) 
4:71 n Jo V2 27r / 

determines (3 at the phase boundary. The resulting curve is shown in the phase diagram 
in Fig. El (label "k = 1"). An expanded plot which reveals more details about the shape of 
the curve is displayed in Fig.lHl Eq. (jSU)) alone with values of (a, z/) that are not restricted 
by Eq. gives the connection between a, v and (3 in the massive phase away from the 
crystal region. 

As in the case /t — > we can perform an analytical near-boundary comparison between 
the crystal and the massive phase. The result is 

*ren,min — ^ren,K=l ~r , , ^ ) 

TT ln(r - To) 

where the relation between ao and v is given by Eq. (jl^ . Since T > Tq, ln(T — Tq) < 
in the crystal phase near the phase boundary, we have again established that the phase 
transition is genuine: In the region where mathematically both phases exist, the crystal 
phase has lower potential. The (T — To)-dependence again shows that the transition is 
second order. 

Limit T ^ 

Many equations simplify remarkably in the limit T — 0, although the ansatz for S{x) is 
basically unaltered. The solution of the T = case has been presented before |22j. Here 
we want to show how to regain the T = results from the general setting. 

We start from Eqs. (jHH) and pn|) where we use asymptotic relations for large z: 

Imlnr(l/2 + iz) ~ z{\n\z\ - 1), Re^(l/2 + iz) ~ In |z| (here, z = P{fx± AA^)/2ti). The 
integrals are best evaluated for fj, > A, where we do not have to account for the absolute 
value in the logarithms. The general result is recovered by taking the (unique) real value 
of the multivalued analytic continuation to fi < A. Standard techniques (e.g. substitute 
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= X and transform the integrals into contour integrals in the complex plane) yield 



Fl,ren(T = 0) 
F2{T = 0) 

^i(T = 0) 



Re In V/i2 - + ^2^2 + J^2 _ ^2 



E 
K 



Re 



A 



1 - 



E 
K 



'A 



-E 



'A 



-— H Re 

2iT An 



/i2 - ^2 + A^K^ - J 11^ - A2 



(52) 



where F/E(x, k) = Jq{1 — t^) ^^"^(1 — n'^t'^)^^^'^ are the incomplete elliptic integrals of the 
first and second kind. For positive A the only solution of -Fi,rcn = F2 = falls into the 
range A'^{1 - k^) < ii'^ < A^ where Eqs. (|S2D simplify to (use EK' + E'K - KK' = 7r/2) 



i"l,ren(T = 0) 

F2(T = 0) 
^i(T = 0) 



ln(A/t) 
E vr 

A^t^ A^ 



All 



27V 



(53) 



Solving Fi ren = -^2 = givcs 



A = - 



(54) 



and the following parameter representation of the grand canonical potential density, 

1 



^ r 



1 

4:71 



27rK2 



2E 

7CK, 



We confirm that A^(l 



K 



< /j!^ < and find that '^rcn < ^E'r 



(55) 

1) and ^'rcn < 



\&ren(^ = 0) = — /X /27r. (For large /x we have ^'rcn = 
consistent with the phase transition at /x = fj,{K = 1) = 2/71. 

Moreover 



-;u2/27r - i2-^/QA7T.) This is 



p(T = 0) 

Fren(T = 0) 



d_ 



1 



2kK 

1 



and 
1 



An 



+ 



nK" 



E 
K 



(56) 



We recover the value of the Fermi momentum pp = at T = from Ref. |22j. 

In Figs, d and |H1 ^ren is plotted together with the self-consistent solution of the homo- 
geneous phases (two curves for the "massive" phase corresponding to the minimum and 
maximum of \E'ren and one curve for the "massless" phase). 

We end this section with a T = analysis of the x-dependent baryon density. From Eq. 
((121) we obtain 



p(a;,T = 0) 



1 



K' 



2kK 2nK 



(57) 
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This expression has the following high and low-density limits: At low density {k — > 1) we 
recover the result for a single baryon 

p(x,T = 0,K^1) ^-^ + T7^, r- (58) 

' ' ' 4cosh2a; 4:C0sh\x/K + K) ^ ' 

At high density {k 0) the total baryon density is constant 

p{x,T^O,K^O)^^. (59) 

Limit a — > 0, perturbative phase boundary and tricritical point 

For all prominent functions in this paper the power series at a = can be given in closed 
form. With c„ defined by 



a;2 



we find Cq — 1, Ci — 1 — C2 = 1 — + 3k,^/S and in general 



E CnC^"'" (60) 



7r/2 n/\/\ / 2\'= 



n 



7T 



We obtain 



^2 = y\ n \ TTT ^V^Re^ 2n, - + — 62) 



n=l 



An^ I (2n - l)!(2n+ 1) V ' 2 2ttJ 



^ 6 2 47r2j (2n + l)! 



Re^(2n,i + ^)-5„,oln£ 



In the following we want to study the tricritical point by expanding the thermodynamical 
potential density to order a^. The result from minimizing \E'rcn with respect to a and k is 
encoded in the fourth order expansions of Fi ren = ^2 = 0. From F2 = we find 



10C(,.') '^^'^-^ + fe' (63) 



with 

G,{k') = )' \r E? (64) 



2^ _ (1 - K^)co - |ci 

(1-«2)C1-|C2 



and the limits 

Gi(0) = | Gi(l) = ^. (65) 
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If the right hand side of Eq. ()63|) is positive, there exists a non-trivial solution for a. The 
equation Fi = determines j3 as function of k and v to 



In 



47r 



1 iz/ 







ReV^(2,i + f; 



with 



and the limits 



f C1C2 



(1 



G2(0) 



Q 2 
The potential density as function of k and v is given by 



(66) 



(67) 



(6^ 



TT 



'I 



1(1 + 



:i - «:2)ci - 



kC2 



ReV'(4 



I 1^ ■> 

' 2 ~^ 27r> 



(69) 

The two phase boundaries of the crystal phase correspond to k = 0, 1. By using Eqs. (jHEI) 
and (jU7|) . one can plot T = l/[3 against jj, = u j [3 ioi u > Vt and construct these phase 
boundaries near the Lifshitz point. The third, "old" phase boundary follows similarly 
from Eq. ()45|) for v < Ut. Finally, in the vicinity of the Lifshitz point, we can expand 
(T, /x) around (Tt,fit) and derive the (approximate) phase boundaries in closed analytical 
form. Starting from Eq. (j66p . we find to second order 



In 



I 



with 



A, 



-AAu - ut) - A2{y - ut) 



(70) 



27r 

5 



— Imip 1, — I 



2 27r 
Im^(3,i + f 



24vr2 Re^(4,i + f^) 



(71) 



The equation for lines of constant n then reads 



T-T, 



AiTt 



UtAi + Tt 



{fi - fit) - 



{A\ - 2A,)Tl 
2{iitAi + Ttf 



(72) 



The linear term is independent of k, so that all curves enter the Lifshitz point with the 
same slope. The quadratic term describes the ^-dependence of the curvature and changes 
sign as a function of k. The limiting cases /t = 0, 1 yield the approximate phase boundaries 
in closed form. Furthermore, Eq. (jUHjl shows that a? is linear in v — z/^. This translates 
into the behavior 

A~(/i-/iO^/' (73) 
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along the phase boundaries and the corresponding critical exponent 1/2 as expected by 
mean field theory. 

Finally, we derive a Ginzburg-Landau type effective action near the Lifshitz point from 
the Taylor expansion of the grand canonical potential in powers of k. For this purpose, 
we compute spatial averages of powers of S and its derivatives (keeping all even terms 
up to order A^). Upon comparing the (analytical) results of such a calculation with Eqs. 
(jUT|) . we find the simple relations 

E 3 

C2 - — Cl 



K 8/14 
E 5 



C3 - T7C2 



{{S') + {{S'f)] 
{{S') + \{{S"f) + h{S\Sr)) (74) 



K ' 16^6 

This enables us to write down a Ginzburg-Landau effective action as follows 



TT^o jj? 1 ^, 

^eff = --T^-^ + —S 



^2 



K4.r) + Re..(l + ^) 



^(S' + (Sr)Re^(2.1 + ^) (75) 



26 

( + -{Sy + 5S\S')A Re^(A,- + 
2ii37r5T4 V 2^ ' ^ ' J ^ V 2 27iT J 

(This result can be applied to the translationally symmetric solution by simply dropping 
all the derivatives of 5*.) At the tricritical point, both the 5*^ and 5"^ + (S")^ terms vanish, 
see Eqs. fl^ and (jlTj) . This is the reason why we had to include terms up to the order 
S^. 

For completeness we report that we have found another self-consistent solution of the Dirac 
Hartree-Fock equation which corresponds to the j = 2 (double gap) Lame equation. The 
scalar density for this solution is 

= 6A^ sn(Ax,.)cn(Ax,.)dn(Ax,.) ^^^^ 

1 + + V 1 — + — Sk'^sii'^^Ax, k) 

where again A and n have to take specific values to ensure self-consistency for given /i 
and T. The solution was found to be analytically accessible, being just slightly more 
involved than the j = 1 case. However, the result showed that (as physically expected 
from the gap-structure) the solution always lies higher in grand potential than the j = 1 
case. Moreover, in general it does not correspond to a local minimum in the space of 
scalar potentials, having unstable directions. Therefore we refrain from writing down 
explicit equations. Instead we show in Fig. |H1 a closeup of the T = plot (Fig. Ej) of the 
grand canonical potential density with the j = 2 case included (in its range of existence 
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5 Relation to condensed matter physics 



In a previous paper, we have pointed out that GN type models can be reinterpreted as 
relativistic superconductors |27j. This is due to a 2-dimensional remnant of the Pauh- 
Giirsey symmetry of massless fermions |2Hll2n|- It allows us to define "particle" and "anti- 
particle" independently for left-handed and right-handed quarks. Models with fermion- 
antifermion pairing (chiral condensate) and fermion-fermion pairing (Cooper pairs) can 
then be mapped onto each other by a "duality" transformation. The phase diagram which 
we have discussed above is equivalent to the phase diagram of a theory with Lagrangian 

C = ^«i#« + ^ + 4V^y , (77) 

provided we replace the chemical potential fi = fiji + fi^ hj the axial chemical potential 
A*5 = fJ^R — fJ^L (see Ref. |2Z|)- The kink-antikink phase of the GN model can then be iden- 
tified with the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase [SHI EI of the dual model. 
Such inhomogeneous superconductors have recently attracted considerable attention in 
the context of QCD (for a current review article, see Ref. [S2])- 

Actually, there is yet another kind of connection between the GN model and non-relativis- 
tic condensed matter systems which is even more surprising. Although the GN model is 
mentioned occasionally in the corresponding condensed matter literature (see e.g. 
l3^). it seems that this close relationship has never been exploited in a systematic way. 
The physics of the problems we have in mind is totally different from the GN model, 
yet the mathematical analogies are striking. Let us discuss a few pertinent references 
which we came across when looking for parallels between the phase diagram of the GN 
model and that of quasi-one-dimensional solid state systems. This may be regarded as a 
continuation of the work of Jackiw and Schrieffer on kinks in condensed matter and 
relativistic field theories to finite kink densities. 

The first example is a 1981 paper by Mertsching and Fischbeck on the Peierls-Frohlich 
model [Sn! • These authors address the quasi-one-dimensional Frohlich model with a nearly 
half-filled band, an electron-phonon system. They are interested in the phase diagram, no- 
tably the transition between commensurate-incommensurate charge density waves, using 
a mean field approximation. There seems to be a mathematical one-to-one correspondence 
between this system and the GN model. The authors of Ref. [32] have also found the 
analytic solution to the mean field equation, guided by the Landau expansion around the 
triple point (which is called Leung point in this context). Crucial for the close corre- 
spondence with a relativistic field theory are evidently a continuum approximation (the 
lattice constant acts as inverse cutoff and is taken to 0) and a linearization of the electron 
dispersion relation near the Fermi surface (simulating "ultrarelativistic kinematics"). In 
appropriate units, the results are identical to ours, as far as we can tell. 

In a subsequent paper, Machida and Nakanishi [SHI used the phase diagram of Ref. |^ in 
a different physics context: They studied the interplay of superconductivity and ferromag- 
netism in ErRh4B4 (Erbium- Rhodium-Boride). Thanks to a number of approximations 
(in particular the one-dimensional band model with linear dispersion) they managed to 
reduce this problem mathematically to the Peierls-Frohlich model. For real order param- 
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eter, their results are again fully equivalent to ours for the GN model, except that now 
one has to use another dictionary: The Dirac equation corresponds to the Bogoliubov- 
deGennes equation, particle/antiparticle degrees of freedom to spin (which exists in a 
quasi-one-dimensional world), chemical potential to magnetic field, baryon density to 
spin polarization. Our three phases (massive, crystal and massless) correspond to their 
BCS, "sn" and normal phases, respectively. The order parameter at T = looks different 
from ours at first sight. However, the two expressions can be converted into each other 
by Landen's transformation for Jacobi functions in the form 



Not only the phase boundaries, but all observables can be identified if one keeps in mind 
the above mentioned dictionary. Nevertheless, there are subtle differences which show that 
one should not take over results blindly: The authors of Ref. fHEj also discuss complex 
order parameters, in particular a helical phase. In the discrete chiral GN model, the 
order parameter is always real. In the continuous chiral NJL2 model, it can be complex 
and we have indeed found a helical phase |,15j, but the results for the phase diagram 
are different. Here apparently the quantitative correspondence between relativistic field 
theory and solid state physics ends due to some differences in dynamics. 

As a third example, we would like to mention the more recent work of Buzdin and 
Kachkachi [3Hj- They derive the Ginzburg-Landau theory for nonuniform LOFF super- 
conductors near the tricritical point in the (T, if)-phase diagram in one, two and three 
dimensions. If we take their result for one dimension and specialize it to a real order pa- 
rameter, we find perfect agreement between our Eq. (fTSj) and their Eq. (3) in appropriate 
units (a discrepancy in the sign of the term must be due to a misprint). We have to 
identify their magnetic field Tio with our chemical potential fi. One can find many more 
similar studies in the literature, see e.g. the recent review article on stripe phases |3CT| . 

Our fourth example is related to solitons in trans-polyacetylene, starting with the seminal 
paper by Su, Schrieffer and Heeger ||41j. In this case there have been efforts to bridge 
the gap between particle and sohd state physics by Chodos and Minakata US] • Since 
these authors start from the old phase diagram of the Gross-Neveu model, they notice 
that there are some discrepancies with the soliton lattice theory of Horovitz jUj. How- 
ever, the analytical results of Horovitz agree perfectly with our kink-antikink crystal so 
that everything falls into place. Let us also mention in passing the results of Okuno 
and Onodera on the coexistence of a soliton and a polaron in trans-polyacetylene [IH] 
which are apparently closely related to the recent work of Feinberg on the analogous 
phenomenon in the Gross-Neveu model. 

Finally, when browsing through the literature on superconductivity, we were amused to 
discover that the "old" phase diagram of the GN model (assuming unbroken translational 
invariance) is ubiquitous in textbooks and articles on BCS theory (see e.g. Ref. (47|j, Fig. 
6.2). It seems to go back to a 1963 paper |48j which predates the original GN paper 
by more than one decade and the first published GN phase diagram by more than two 
decades. 





(78) 
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Summarizing, we find it gratifying tliat tlie GN model witli its simple Lagrangian gives 
rise to such a physically rich phase diagram. Its structure seems to be mirrored by a 
variety of quasi-one- dimensional condensed matter systems. With hindsight, one might 
ask the question why it took much longer to solve the same mathematical problem in 
relativistic quantum field theory than in condensed matter physics. We feel that this is 
due to the almost exclusive use of path integral methods in particle physics. Apparently, 
these methods are not yet sufficiently developed to yield non-trivial saddle points with 
unexpected symmetry breakdown as in the present case. By contrast, canonical mean- 
field calculations offer a rich spectrum of approximation methods (variational, numerical 
etc.) which can be successively refined and used to nail down eventually the exact solution 
as in the case at hand. It is probably not accidental that we arrived at the phase diagram 
of the GN model using Hartree-Fock methods akin to the standard tools of condensed 
matter physics. 
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upper band 



lower band 



2K P 2K 



Figure 1: The dispersion relation pij) for k = 0.8. The Bloch momentum p is projected 
onto the Brillouin zone ranging from — 7r/2K to +7r/2K. The bold line marks the branch 
for which the slope of the dispersion relation is positive. Notice the symmetry p —>■ —p. 
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0.6 0.8 1 1.2 1.4 1.6 

Figure 4: Contour lines k — const, in the crystal phase (closeup). The phase boundaries 
are k, — and k, — 1, adjacent lines differ by Ak, — 0.05. All lines end at the tricritical 
point. At T = we have Ak, — 1 and fj, — 
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Figure 5: Phase diagram of the GN model, computed with the help of Eqs. (jlH), 
(Uni); (j^nj) . All phase boundaries correspond to second order transitions, Pl is the Lifshitz 
point determined by Eqs. (j^5|) . ()47|1 . 
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Figure 7: The grand canonical potential density ^ for self-consistent solutions of the 
Dirac-Hartree-Fock equation at T = 0: (a) masslcss Fermi gas (m = 0), (b) and (c) the 
minimum and maximum solution for the massive Fermi gas (m = 1) and (d) the crystal 
solution with the Lame potential with single gap (j = 1). Notice that the crystal lies 
below the massless Fermi gas. At = 2/7r, a second order phase transition between the 
massive Fermi gas and the crystal occurs. 
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Figure 8: A closeup of Fig. [7| showing \1/ at T = for all known self-consistent Hartree- 
Fock solutions. Here the (unstable) double gap Lame solution (j = 2) has also been 
included (e). 
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